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Summary 

The Euler equations of gas dynamics have some very interesting properties in that the flux vector 
is a homogeneous function of the unknowns and the equations can be cast in symmetric hyperbolic 
form and satisfy the entropy conservation. Since the Euler equations are the moments of the Boltzmann 
equation of the kinetic theory of gases when the velocity distribution function is a Maxwellian, it would be 
interesting to look for the relation between the symmetrizability and the Maxwellian velocity distribution. 
The present paper precisely shows this relationship. The functions that symmetrize the Euler equations 
are density and mass flux, which are integrals of the Maxwellian distribution multiplied by unity and 
molecular velocity, respectively. The field vector and the flux vector in the Euler equations are then the 
gradients of these functions with respect to new transformed variables q. In terms of these transformed 
variables, the Euler equations assume the symmetric hyperbolic form. 

The entropy conservation is in terms of the //-function, which is a slight modification of the //-function 
first introduced by Boltzmann in his famous //-theorem. The modified //-function is equal to the negative 
of specific thermodynamic entropy of the gas. Further, it is an integral of the Legendre transform of the 
Maxwellian distribution with respect to q, thus establishing its convexity. It is therefore possible to 
design a numerical method using the convexity of H so that the total H in the computational domain 
will monotonically decrease with time in conformity with the Boltzmann //-theorem. In view of the 
//-theorem it is suggested that the development of total //-diminishing (THD) numerical methods may 
be more profitable than the usual total variation diminishing (TVD) methods for obtaining “wiggle-free” 
solutions. 

Introduction 

The Euler equations of gas dynamics can be obtained as the moments of the Boltzmann equation of 
the kinetic theory of gases (ref. 1) provided the velocity distribution function, which is the basic unknown, 
is a Maxwellian distribution. The collision term in the Boltzmann equation vanishes when the velocity 
distribution function is a Maxwellian. The Euler equations therefore are the moments (the moment 
functions being collisional invariants) of the collisionless Boltzmann equation that are linear hyperbolic 
equations with a wave velocity independent of space and time. Because the wave velocity does not depend 
on space and time, the collisionless Boltzmann equation can easily be cast in the strong conservation law 
form. Taking the moments of this equation therefore yields the Euler equations in the strong conservation 
law form. This fundamental connection between the linear hyperbolic Boltzmann equation and the Euler 
equations has been used in reference 2 to construct a new class of upwind methods for the numerical 
solution of the Euler equations. These methods, called “kinetic numerical methods,” are based on the 
principle that moments of every upwind method for the Boltzmann equation yield an upwind method for 
the Euler equations. Central to the moment-method strategy is the Maxwellian distribution. 

The Euler equations can be cast in two forms: the strong conservation law form and the symmetric 
hyperbolic form. The strong conservation law form, which is a consequence of the conservation of 
mass, momentum, and energy, is the basis for constructing fully conservative methods. However, the 
wave nature of the Euler equations is not apparent when the equations are cast in this form. The 
symmetric hyperbolic form studied by Harten (ref. 3) makes the hyperbolicity of the Euler equations 
very transparent. This effect is due to the symmetry of the Jacobian matrices occurring in the symmetric 
hyperbolic form. The symmetrizability is a very important property of the Euler equations. The excellent 
report by Harten (ref. 3) shows the connection between the symmetrizability, entropy function, and Roe 
linearization. According to the Godunov theorem (ref. 3), symmetrizability implies the existence of the 
entropy function which, in turn, according to the Harten-Lax theorem (ref. 3), implies that the equations 
admit Roe linearization. Hence, the symmetrizability of the Euler equations allows the equations to 
be locally linearized so as to preserve hyperbolicity and conservation. Symmetrizability is thus a very 
important property of the Euler equations both from theoretical as well as numerical points of view. In 
fact, Abarbanel and Gottlieb (ref. 4) have used the symmetrizability property to analyze rigorously the 
splitting algorithms for the Navier-Stokes equations. Further, it is possible to construct upwind methods 


(ref. 5) based on Roe linearization, for example, that lead to the improved structure of iteration matrices, 
i.e., diagonal dominance. 

With the symmetrizability and the entropy condition being very important properties of the Euler 
equations, it would be very interesting to seek their basis in the fact that the Euler equations are 
the moments of the Boltzmann equation. The motivation of the present paper lies precisely in this 
observation. The present paper shows that the Maxwellian distribution plays a very important role in 
the theory of the Euler equations. The integrals of the products of the Maxwellian distribution and first 
two collisional invariants (corresponding to mass and momentum conservation) are the functions that 
accomplish symmetrization. The entropy condition is related to the Boltzmann //-theorem. This relation 
is to be expected especially when it is noted that according to the analysis of the kinetic numerical method 
of reference 2, the collision phase decreases //, but the convection phase conserves H. 

It is believed that the results of the present paper will be valuable in constructing a new class of upwind 
methods using the moment-method strategy for obtaining “wiggle-free” solutions. For the Euler equations 
the moment-method strategy is based on the Maxwellian distribution and the Boltzmann equation. 
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Jacobian matrix of flux g with respect to w 
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constants dependent on 7 

specific total internal energy per unit mass 

Maxwellian distribution 

contracted Maxwellian distribution defined by 
equation (24) 

arbitrary distribution 

flux vector in Euler equations 
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degrees of freedom 


2 


M 0 

M+,M~ 

M 0 
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Mi 

dM 0 ,dM i 
dM x , dM y 

P 
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integral of F with respect to molecular velocity 
and I 
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v being half-intervals (0,oc) and (— oo,0) 

any one of M 0 , M+, and M~ 

integrals of v\F, and v$ F, respectively, with 
respect to v\, V2 , 1^3, and I 

integral of vF with respect to v and / 

differentials of M 0 and M \ , respectively 

differentials of M x and M y , respectively 

matrix appearing in symmetric hyperbolic form 

pressure 

quadratic form 

transformed variables in symmetric hyperbolic form 

element of q 

value of q at time level n 

values of q at neighboring points 

gas constant per unit mass 

thermodynamic entropy per unit volume 

temperature 

time 

fluid velocity 

molecular velocity 

molecular velocity in z-direction 

field vector in Euler equations 

elements i and j of field vector w 

field vectors with u-integration over half-intervals 
(0, 00) and (—00, 0) 

differentials in equations ( 53 ) 

coordinates 

square of thermal speed 1/(2 RT) 
ratio of specific heats 
general vector 
transpose of vector rj 
coefficient between 0 and 1 
mass density 


* 


collisional invariant vector 


vJ/ ? ith component of vector *P 

Abbreviations: 

THD total //-diminishing 

TVD total variation diminishing 


Homogeneity and the Maxwellian Distribution 


The one-dimensional 
vation law form 


unsteady Euler equations of gas dynamics can be written in the vector conser- 


gw dg 
dt dx 


= 0 


(i) 


where 



" p " 


pu 


pu 

w = 

pu 
. pe _ 

g = 

p + pu 2 
_{pe + p)u _ 

— 

p + pu 2 

.^T pu + ^pu 3 _ 


( 2 ) 


Here, p, u , and p are, respectively, the mass density, fluid velocity, and pressure, and e is the specific total 
internal energy given by 


pe = 



1 2 

+ 2^ 


(3) 


Equations (1) can be obtained from the Boltzmann equation of the kinetic theory of gases (refs. 1 and 2) 
and can be written as 


T dF dF\ ^ 

'"■ar + '’W =0 


(4) 


where VP denotes the moment functions 1, and I + (v 2 /2) corresponding to the collisional invariants, 
and F is the Maxwellian velocity distribution 


F 



_ „)2 _ _ 


(5) 


Here, (3 = l/(2/?T), R is the gas constant per unit mass, I Q = (3 — 7)/ [4(7 — 1)/?] that denotes the 
internal energy due to the nontranslational degrees of freedom, v is the molecular velocity, / is the 
molecular internal energy variable, and the inner product is given as 

/ oo roo 

dv dI{*F) (6) 

-00 Jo 


The formula for the Maxwellian distribution (eq. (5)) contains two independent variables, the molecular 
velocity v and the internal energy variable I . The latter variable is required to ensure the existence of 
additional degrees of freedom that are necessary to satisfy the constraint 


I+-z,F) = 


2(7 


1 ^ + 2 ^ 


The parameter 7 is assumed to be a constant for the analysis in the present paper. Notice that 7 0 , which is 
the average energy in nontranslational degrees of freedom, will be positive if 1 < 7 < 3. For a monatomic 
or polyatomic perfect gas, 1 < 7 < 5/3 and, hence, I Q will always be positive for such gases. 
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The unknown field vector w and the flux vector g are related to F by the equations 


w 


= (9<F) = I 


1 

V 

/+£ 


Fdvdl 


2 J 


g = (v¥,F) = J 


1 

v 


Fdvdl 


(7) 


( 8 ) 


. / + \ . 

where only one sign of integration is displayed for brevity. The homogeneity property of the flux vector is 

Aw = g (9) 


vheve A — <9g/<9w. In terms of F we have 


( Aw) f = £ A ijm , = E ||»> = / »*, fe 

3 3 J \ 3 , 


dv dl 


Hence, the homogeneity property (eq. (9)) follows if we can show that 

dF 


( 10 ) 


By using the chain rule of partial differentiation, we get 

dp \ dF 

' w >d^r~d-p \^ Wj 

3 


Using 


gives 


dF dF ^ dp \ dF du \ dF ^ df3 

? ■ = Tp l j + Tu (E I + M I 


80 


3 dwj 


w 2 1 of 01^3 w 2 

P = wi u =— * = 2(Tf-l)| — 


dp 


= 1 


dw\ 

du u 

dw\ p 


*2 2 


90 0 

= - - n- U 


p p 

From these equations it follows that 


E 


w,- 


dp 

2 dwi 


u> l 0 

du 1 
dw2 p 

d0 _ 20 2 (^-1 )u 
dw2 p 




dp_ 

dw$ 

du 


= 0 


= 0 


re 


dw 3 

dl _ 2^ 2 (7-l) 

dm3 p 

d0 


J dw~ 


= 0 


( 11 ) 


and, hence, equation (11) yields the desired equation 

dF 1 dF 

1 
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It is interesting to note that the homogeneity property (eq. (9)) depends on the defining relations (eqs. (7) 
and (8)) and on equation (10). Hence, the velocity distribution need not be a Maxwellian for the validity 
of equation (9); it is enough if equation (10) is true. 


Symmetrization 

The equations 


dw . dw 

~di +A l^~° 


are said to admit a symmetric hyperbolic form if they can be transformed to 

-s-a- 


( 12 ) 


where P and B are symmetric matrices and P is positive definite. We will accomplish the symmetrization 
of the Euler equations by obtaining the scalar functions M 0 and Mi such that 


dMn 

dM 0 

dM 0 

(13) 

CO 

II 

3 

™2 = ~ 
dq2 

" 3 = a®, 

dMi 

91 ~ a qi 

dM\ 
92 dq 2 

SM, 

93 = a® 

(14) 


where qi , q 2 , and 93 are yet to be determined functions of p, u, and (3. Assuming that M 0 and Mi have 
been obtained, the time derivatives of w and the space derivatives of g are given by 


dwj 

dt 

dgj 

dx 


dx dq 

Hence, the Euler equations (eq. (1)) transform to 

d 2 M 0 dqj 


d dMp _ dHfo_dqj_ 
dt dq{ — dqi dqj dt 

d dMi _ ^ 3 2 Mi dqj 


E 


T* dq t dqj dt dqj dqj dx 


— ' dqj dqj dx 


9 2 M\ dqj _ ^ 


(15) 


3 J 3 

and the required matrices P and B are then given by 

P = 


' d 2 M 0 ‘ 

O 

' <9 2 Mi ’ 

dqi dqj 

15 — 

% dqj 


(16) 


We will now prove that the desired scalar functions M 0 and M\ accomplishing symmetrization are, 
respectively, 


M 0 = 

roc roc 

/ dv dl (F) = 

[ Fdvdl 

(17) 


J —00 J 0 

J 



r 00 ro 0 

r 

(18) 

Ml =J 

I dv dl {vF) = 

—00 Jo 

/ vFdvdl 


The differentials dM 0 and dM\ are then given by 


dM 0 

_ f 

[11 

_dM\ 

-J 

£>' 

1 


dF dv dl 


(19) 
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Using 


dF 

dp 


F 

P 


^ = 2 :«» - u)F 


dF 

d/3 


A (r u) 2 4(7-1) 

L 2/3 (U U) 3 - 7 J 


( 20 ) 


we obtain 


— = 1 - 2(3(v - u)du + 


dF _ dp 

T~~ P 

dp 

P 


A _ {v _ u) 2 _ linilj 

L 2/3 [v u > 3-7 


d/3 


KlU O O 

= — - 2/3w du -u z d(3 + —d/3 + (2/3i> du + 2 uu d/3) - 


' 4(7-1) 

. 3-7 


7 + v 2 


d/3 


Using total differentials we obtain 


dF 

F 


d ^ln p + | In (3 - (3v?^J + v d(2/3w) 



d( — 2/3) + — ^ 37) 7d/3 

3-7 


( 21 ) 


where the identity 


4(7- 1) 2 2(5 -37) 

3-7 3-7 


has been made use of. Substituting for dF in equation (19) gives 


dM 0 
dM i. 




^ ln/3 — /3u 2 


^ + vd{2f3u) 


+ 


M)«- 


2/3)+^-Midf3 
3-7 


Fdvdl 


( 22 ) 


For molecules with additional degrees of freedom, i.e., 7 ± 5/3, more manipulation is required to absorb 
the 7 d/3 term with other terms of equation ( 22 ) in order to express the integrands as sums of products of 
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collisional invariants and perfect differentials. This manipulation is accomplished by using the following 
equalities: 


2(5 - 3 7 ) 


( Id0)FdvdI 


2(5 - 37 ) 


f Io(~Fd(3)> 


2(5-37) 3-7 

3-7 4(7-1)/? 


5 — 37 {dp 


{F dp) dv 


f 5-37 
J 2(7-1 


F dvdl 


where F is the contracted Maxwellian given by 


= J F dl = p\j^exp [— P(v — u) 2 


By following a similar procedure we can easily prove 


m 


2(5 - 37 ) 
3-7 


'vIdP(F)} dvdl = [ v^- FdvdI 

2 n - i) J P 


Using equations (23) and (25), equation (22) reduces to 


J l^j ^lnp + - In/? - j + vd(2Pu) 

+ { I+ j) d{ ~ w) + ^=T) d{ln/3) FdvdI 
J * d (in p + ^ - n ^ - /?u 2 ) + v d(2Pu) 


+ I 7+ — 1 d{-20) FdvdI 


The transformed variables q \ , < 72 , and <73 can therefore be defined as 

In 3 

qi = In p H - Pu 2 q 2 = 2/3u q 3 = -2 P 

7-1 

Equation (26) then assumes the simple form 


11 v z \ 

v dq i + v dq 2 + I / + — 1 dq 3 F dv dl 


Equation (28) then implies 


^ dM ° = [ 

dq t J 

<9 Mi f 

~ dq t 


Vi FdvdI 


v^i FdvdI 


Equation (28) is clearly the crucial relation in accomplishing symmetrization for it reveals that the 
gradients of scalar functions M a and M\ are, respectively, equal to w and g. It is therefore worthwhile to 
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see more transparently the connection between equation (28) and the Maxwellian distribution. In terms 
of 9 i> 92> and 33 defined by equation (27), the expression for dF reduces to 


dF , , / u 2 \ J 

— =dqi+vdq 2 + [I + — ] dq 3 + 


= dqi + v dq 2 + I I + — dq$ + 


= dqi + v dq 2 + 


K) 


2(5 - 3 7 ) 
3-7 

5-37 
4(7 - 1)/? 


I dp- 


5 — 37 dft 


dqz - 


5-37 

3-7 


Idq-i 


dqs + ^-^ilo-Vdqs 
3-7 


The validity of equation (28) is obvious when examining equation (31). Thus, 

ro o 

/ (7 0 — 1)F dl = 0 

^0 


(31) 


(32) 


The above analysis makes use of the Maxwellian distribution to establish the validity of equations (29) 
and (30). Once qi, q 2 . and <73 are defined by equation (27), the validity of equations (29) and (30) can be 
directly verified from the definitions of M 0 and M\. Equations (17) and (18) give, respectively, 

M 0 = p M\ = pu 

From equation (27) we obtain the total differentials 


j _ dp dp 

dq ' = J + h ^ 

dqi = 2u d/3 + 2ft du 
dq 3 — —2 dft 

The expressions for dq 2 and dq 3 yield 


— u 2 dft — 2 ftu du 


du — 


dq 2 + u dq 3 

2ft 


Substitution of du and dft into the expression for dq\ gives 

dq ' = 7 + (-t) + “ 2 (tO “ + udq3) 


dp 

P 


1 


iP 

2 + 2(7-1 )p 


dq 3 — u dqi 


Hence, 


dM 0 — dp ~ pdqi -b pu dqi + 


pu * 


+ 


2 ' 2(7 - l)ft 

which implies that w 1 = dM 0 /dq Similarly, for dM\ we get 

dqi + 


dqz 


dM\ = pdu ■+- u dp = p 


udq 3 J 
— +u \ 


dMi = pu dqi + 


2/3 

+ pu 2 ) dq 2 + 


(v> +pu2 ) 


pdqi + pudq 2 + 
3 


2(7 - 1 )P 


+ 


pu 


dqz 


~fpu 


+ 


pu" 


2(7-1)/? 2 


dq 3 


rfAfi = gi d</i + 32 <?32 + 33 <% 


(33) 


(34) 


(35) 


(36) 
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Thus far we have shown that the derivatives of the scalar functions M 0 and M\ with respect to q{ are, 
respectively, equal to and g±. Performing differentiation once again yields the elements of the matrices 
P and B. The symmetric hyperbolic form for the Euler equations can then be written in the expanded 
form: 


P pu 


a 

'In p+£4-/J« a - 

pu p + pu 2 

o P u+e£ 

0 

dt 

2/3 u 

+ apu+ e£ 

c^+opu* + £|1_ 


-2/3 


pu 

p + pu 2 

cpu + e^- 


’In P+^-0u*- 

p + pu 2 

3 pu + pu 3 

-£ + ( CT+ 3 )pu 2 + £ «l 

d 

dx 

2/3 u 

pu 3 

_ crpu + 

£h! + ( <7+ 3 )pu 2 + £ «1 

(<t 2 + a) 2^ + (a + 1) pu 3 + _ 


-2/3 


I 

i 

I 

! 

| 

(37) ! 


where a — 7/(7 — 1). It is interesting to note that the scalar functions M 0 and M\ symmetrizing the 
Euler equations are, respectively, mass density p and mass flux pu. The Euler equations (eq. (1)) are in 
conservation form (w-representation) and transform to symmetric hyperbolic form (q- representation) in 
terms of new variables q. The w-q transformation is given by 


dM 0 _ dp 
3q <9q 


3 Mi 
g “ dq 



The w-representation reflects the physical principle of conservation, whereas the q- rep resent at ion clearly 
reflects the hyperbolicity. The mass density and the mass flux are at the root of the above transformation. 
This preferential role played by the mass density and mass flux is physically due to the fact that mass is 
the carrier of the momentum and energy. 

Finally, we show that equations (13) and (14) are valid even when the integration with respect to v in 
equations (17) and (18) is over either of the two half-intervals (0, 00 ) and (— 00 , 0). The integrations over 
half-intervals are required when the Maxwellian distribution is split into two truncated Maxwellians — one 
corresponding to particles with positive v and the other corresponding to negative v. Such a splitting 
gives rise to splitting of the flux vector g into g + and g _ . The split Euler equations (corresponding to 
the splitting of g) can also be cast in the symmetric hyperbolic form using the integrals of the truncated 
Maxwellian. To demonstrate this property, we offer the following definitions: 


r w+ n 

w t 


*1 

l 

roc roc 

= dv dl 

v 2 

Jo Jo 

[i+ v 4\ 


91 

92 


r 00 roc 

1 

— 1 dv dlv 

V 0 

Jo Jo 

1 

+ 


'm+' 

f OO r 00 I 

r 

M+ 

^0 

0 

11 

V 


Proceeding as before we obtain 


~dM + 1 f°° f°° , 

r 

r 

Wt\ -Jo * 7 . dI 

_v_ 

F 


dqi + v dq 2 + I I + y 1 dq% 


+ ^ 1 (Io-I)dq 3 
3-7 


(38) 


(39) 


( 40 ) 
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In view of equation (32), the above equation simplifies to 


dM+' 

dMf 




dq\ + v dq 2 + 



dqz 


Equation (41) then immediately implies that 



dM+ 

dq t 



<9M+ 

dqi 


(41) 


(42) 


Thus, when w-integration is over the half-interval, the transformed variables , < 72 , and <73 are the same 
as before but the scalar functions are M+ and instead of M a and Mi . 

By proceeding on similar lines we can establish the following results. Let 



(43) 


(44) 


(45) 


The integration with respect to v in equations (43), (44), and (45) is over the half-interval (— oo,0). Then, 
the following relations hold: 


_ dM~ 


w i = 


dqi 


_ dM7 

rt L 


Si = 


dqi 


(i = 1,2,3) 


(46) 


Here < 71 , < 72 , and 93 are the same as before and are defined by equation (27). 
To summarize, the important result of this section is that 


/ dF dF\ v 2 


can be cast in the symmetric hyperbolic form (eq. (15)) irrespective of whether the integration in the 
inner product with respect to v is over the doubly infinite interval (— 00 , 00 ) or the half-intervals ( 0 , 00 ) 
and (—oo,0). The transformed variables appearing in the symmetric hyperbolic form are qq, < 72 , and <73 
defined by equation (27), and the scalar functions accomplishing the symmetrization are mass and mass 
flux. 


Positive Definiteness of P 

The scalar function M 0 is a convex function of < 71 , <72 > an d 93 if the matrix 


P = 


d 2 M Q 

dqi dqj 


11 



is positive definite. The positivity property of P can be easily proved when the integration with respect 
to v is over the full space (— 00 , 00 ). We first obtain expressions for d 2 M 0 jdq{dqj and then determine 
whether the determinants corresponding to the leading minors are positive. For this purpose we have 


d 



— dp — p dq\ + pu dq2 + 


pu 2 t p 

T + 2(7 - 1)0 


dq:i 


(47) 


d (^t) =dM = dMl 

= pudq l + (^~+ pu 2 ^j dq 2 + 


Ypu pvr 

2(7 - 1)0 + 


dq 3 


'(S- 


pu * 


+ 


2 1 2(7 -l)/3 
2 1 


+ 


u 

Y + 2(7 -1)0 

U 2 1 

T + 2(7 -1)/? 
P dq 3 


dp — 


L2(7-1)/? 2 J 


dP + pu rfw 
2 


+ /m 


pd<?l 4- pudq2 4- 
dq2 + u dqs 


pu * 


4- 


2(7 -1)0 


dq 3 


2(7-l)0 2 2 20 

where use has been made of equations (33) to (36). Rearrangement of terms gives 




pu 


+ 


4- 


2 ’ 2 ( 7 - 1 )/? 
7 P 




2 4 

7 /nr ( /nr 
(7 — l ) 2 402 + 7 - T 20 + ~ 


^3 


The second derivatives of M 0 are then given by 


d 2 M 0 


; 2 ~ P 


d 2 M 0 d 2 M Q 


dqj 


0 2 M o d 2 M 0 pu 2 


dq x dq 2 dq 2 dqi 

P 


= pu 


+ 


dqi dq 3 dq 3 dq Y 2 2(7 - 1)0 

d 2 M Q p 2 d 2 M 0 d 2 M Q 

1 - ni‘ 

w 


03 - = - +p “ 


n 

dHdg _ 

dq 2 (7 - l) 2 40 2 n 7 - 1 20 1 4 


d 2 M 0 d 2 M G ^ 7 pu pu 3 

5^3 d <73 d<7 2 7 - 1 20 + 2 


p pu 4 

*• -< 


The determinants corresponding to the leading minors are given by 

d 2 M 0 


dq 2 


= p > 0 


(48) 


(49) 
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d 2 Mn 

~dvt 


d 2 M„ 

dq l dq2 


det P = 


d 2 M„ d 2 M n 

dq? dqi dqf 


P 

pu 


= P\pu + 

pu 

23 + P u2 
pu 2 i p 7 pu , pu 3 

+ 2(7— 1)3 7-1 2/5 + 2" 


2/3/ 


**-*>» 




2(7— 1)3 7-12, 

A sequence of elementary transformations gives 


2(7-1)/? 

7 QU 

7-123 

+ 7^1^ + 


_X_^u + ^ 


det P — 


P 

0 


pu 


2(7-1)/? 

S 


+ T 


2(7— 1)3 TT-12? 

A 


P PU 

0 2)3 

0 ^ 


/9 2 U 2 


2(7-1)/? 


7 p , 27—1 pu 2 
(7-I) 2 4^ + 7-1 43" 


23" + 4(7- 1)3 2 


+ 


/> 2 U 2 


4 /? 2 8(7 - l)/? 3 6/? 2 

>0 (7 > 1) 


8(7 - 1 )/? 3 

Hence, all the leading minors of P have positive determinants and P is therefore a positive matrix implying 
the convexity of M 0 . It is interesting to note that the only restriction on 7 in the above analysis is that 
7 > 1. 

The proof of positivity of P for the cases when the integration is over half-negative and half-positive 
intervals can be constructed on the above guidelines. The algebra becomes very involved and the main 
thrust of the present paper, that many results become transparent in view of the Maxwellian distribution, 
is then defeated. We will now demonstrate the positivity of P in all the cases using the fact that M 0 is 
an integral of the Maxwellian. To this end, we first write M 0 as 




= f b dv f° 

Ja JO 


dIF 


(50) 


where the limit (a, 6) could be any one of the intervals (—00,00), (0, 00), or (— 00, 0). Using equation (31) 
gives 


ro ro 
0=1 dv I 

Ja Jo 

= ( b dv r 

Ja Jo 


dl ( dF ) 


dIF 


dqi+vdq2+ I + — } dq 3 + 


5 — 
3-7 


(■ Io~I)dq 3 


In view of equation (32), the expression for dM 0 becomes 

rb 


fb f c 

dM 0 = I dv I 
Ja Jo 


dIF 


dqi + v dq 2 + I / + — I dqz 
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Taking the differential once more yields 


O - fb roo 

d z M a = dv dl ( dF ) 

J a Jo 


dqi + v dq 2 + I / + — 1 dq 3 


Substitution of dF from equation (31) into equation (51) further gives 


Pm, 


ro re 
0=1 dv I 
Ja JO 


dl dqi + v dq% + ( / + — I dq 3 


dqi + v dq 2 


+ M + — dq 3 + 


5-37 


3-7 


(■ Io~I)dq 3 


Let 


dX = dqi + v dq 2 + ( I + — ) dq 3 


dY = dqi + v dq 2 + ( / + 


dq3 + — — — (I 0 - I) dqz 


3-7 


= dqi + v dq 2 + 
Equation (52) then reduces to 


2(7-1) P 5 — 3 7 
3-7 + 2 + 3 - 7 0 




d l M, 


rb rc 

0=1 dv 
Ja Jo 


dl dX dY F 


Noting that 


dY = dX+%^l{I 0 - I) dq z 


3-7 


dX = dY — ^ — —{I 0 - I) dq 3 


3-7 


( \I 0 -I)dq 3 dX 


equation (54) yields 

pm 0 = [ b dv r d i(dxfF + [ b dv r d i 

Ja Jo Ja Jo L 3 — 7 

and, alternatively, 

d 2 M Q = [ b dv rdI(dYfF- ("dv f°°dl (/ — I 0 ) dq 3 dY 

Ja J 0 Ja Jo L 3 — 7 


(51 


(52) 


(53a) 


(53b) 


(54) 


Substituting for dX and dY from equations (53a) and (53b), respectively, in the second term of the above 
formulas for d M 0 and using the easily provable relations gives 


roo /* oo 

J o (I-Io)FdI = 0 y o /(/ - Io)FdI = I 2 F = J I 2 F dl 
The two equations for d 2 M Q then simplify to 

d 2 M° = £ dv J™ dl (dX ) 2 F — J b dv j™ dl (dqi) 2 F 


(55a) 



d? M a = f b dv [°°dI{dY) 2 F+ f b dv l°° dl [ 2(7 

/ a io fa fo L (3 - 7) 

The quadratic form Q corresponding to the matrix P where 

n = „t ( d 2 M 0 \ 

can be obtained from equations (55) by replacing dqi by rfo. We then obtain 


Q = f dv f 

Ja JO 

-/>/„' 


dl J)i+Vl) 2 +U + y % 


dl 


5-37 

L 3 - 7 


{lorn? 


/ o roc 

dv j dl Irj i + u?72 + 
Ja Jo 


— 1) r , V 2 , 5-37 , 
3-7 2 3-7 ° 




(3-7) 2 


(^3) 2 


(55b) 


It is obvious that Q > 0 as long as 1 < 7 < 3. Hence, P is a positive matrix. Two expressions for Q were 
derived to show that irrespective of whether 7 > 5/3 or 7 < 5/3, we obtain a positive expression for Q. 
If 5 — 37 < 0, the positivity of P follows from the first of the above two formulas for Q , and if 5 — 37 > 0, 
then the positivity of P follows from the second formula. 

The above proof regarding the convexity of M 0 (or equivalently positive definiteness of P) rests upon 
the equation dF — F dY , which in turn is a consequence of the Maxwellian distribution. Once again the 
connection between the convexity of the function symmetrizing the Euler equations and the Maxwellian 
distribution is obvious. 


Entropy Function and the if -Theorem 

The Boltzmann if-theorem has been described in the kinetic theory of gases as the bridge connecting 
the thermodynamics and the statistical mechanics of particles. Briefly stated, the theorem says that the 
//-function (ref. 1 ) defined by 

H = /(/In /) dv 1 dv 2 dv 3 

monotonically decreases with time as a homogeneous gas in statistical nonequilibrium evolves to equilib- 
rium. In case of spatial inhomogeneity, the theorem states that 


where H v ^ is the flux defined by 


dt ^ dxi ~ 

l 1 


H vi - J Vi(f In /) dv 1 dv 2 dv 3 


(56) 


It is therefore natural to define an entropy function for the Euler equations using the above definition of 
H with the arbitrary distribution / replaced by the Maxwellian. We then obtain H by 


and //-flux by 


H = H{F) = j j (Fin F) dvdl 
H v = J J v(F\nF) dvdl 


(57a) 

(57b) 
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When F is a Maxwellian, In F is given by 


In F — In p + - In 0 + In 


4(7 - 1) 
Lv^(3~7)J 


- 0(v - u ) 2 - 4(7 - ~ } PI 

J ~ 1 


which can be equivalently written as 

2(5 - 3-y) 


In F - - ^ 01 = |ln p + £ In 0 + In 


3-7 


3 

2 


4(7 - 1) 


n/^(3 - 1) 


— (3u 2 


+ (2 (3u)v + (—2/3) I I + — 


( 58 ) 


The right-hand side of equation (58) is a linear combination of collisional invariants and, hence, 
substituting equation (4) into (58) gives 


// 


hr- a 

3-7 


1 / dF dF\ 
. \~dt +V dx) 


dvdl — 0 


(59) 


We will show that equation (59) is the basis of the entropy conservation. Considering first that 7 = 5/3, 
equation (59) simplifies to 

j j' aF {ji +v %) dvi,=o < 6o) 

Hence, in view of equation (60), 

y / (i + 4) (FI " H / / (1 + ln f 1 (i- 

In terms of H and H v , the above equation becomes 

dH d H v 


+ v 


5tF s 
dx . 


dvdl = 0 


dt + dx 


0 


(61) 


which is the entropy conservation for a smooth solution. When 7 ^ 5/3, the gas has additional degrees of 
freedom apart from the translational ones. The In F appearing in the definition of H is then no longer a 
linear combination of collisional invariants and, consequently, equation (60) will not be valid. Therefore, 
a slight modification in the definitions of H and H v is required. To obtain these new expressions for H 
and H v we can proceed in the following manner. Equation (59) can be written as 


II 


1 + hr 

3-7 


dF dF\ , JT n 

— + v —j dvdI = o 


(62) 


Using d(F In F) = (1 + In F) dF gives 

//*"(: I + 4) (F F) - / M a I + vd £) dvd, = 0 (63) 

The “trick” lies in transforming the second term in equation (63) as an integral of a perfect differential. 
We note that 

I/' 57 (f + "f ) -II [!<*"■>■ + 4H dvii 

-/Mf^l )" 7 
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where uq is equal to p, , or w ^ , and g\ is equal to gq, or g ^ , depending on the interval of v 
integration. Further manipulation gives 


/MS- a** 

= 3FFli){ (1 + 1 ” <J) (^ 1+ & )" [l ( * ll " ffl + ^ Slln ' ,) ]} < 65 > 

When the integration is over the full interval, 


dw x dg x dp d 

+ rn + di [,m] = 0 

When the integration is over the half-interval, split Euler equations 

dw± 


+ 


= 0 


dt ' dx 

will have to be imposed. Keeping this fact in mind, equation (65) gives 

II 0I {% + "ID ivd ’ = “4^1) [I** 1 M) + §i [h 


Combining equation (67) with equation (63) gives 

//e *-y [■'*»■» 

The //-function and //-flux can hence be defined, respectively, as 


Fdvdl = 0 


(66) 


(67) 


( 68 ) 


H = 

H v = 


II 

II 


dv dl 


dv dlv 


F,nF+ ^~jj Fi,, ^ +c ‘ F 

FlnF + F\n(3 + C 2 F 

47-1) 


(69) 

(70) 


where C\ and C 2 are constants independent of v, /, p, ft, and u. In view of the definitions given in 
equations (69) and (70), equation (68) can be cast in the compressed form 


dH dH v n 

1 - = 0 

dt dx 

which is the entropy conservation. The thermodynamic entropy per unit mass for a perfect gas is 


S = 



In/? 

+ Constant 

7-1 


The appearance of a constant in the above formula is due to the fact that the entropy is indeterminate 
within a constant. The //-function defined by equation (69) after integration with respect to v and I 
becomes 


H = 



In ft 
7 - 1 


+ C\ — - + In 


' 4 ( 7 - 1 ) 1 1 

.v^(3-7)J/ 
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which is therefore negative of thermodynamic entropy per unit volume. Since the //-function is a measure 
of the information content of a distribution, it will be negative of entropy. The modified //-function is 
thus a physically meaningful quantity. 

Because of the appearance of an additional term involving Find in equation (69), the convexity of H 
is not obvious. We will now show that H is a convex function of < 71 , < 72 , and < 73 . The proof rests upon 
showing that H is the Legendre transform of M 0 ; that is, 


TT dM a , dM 0 , dM a 
H = qi — — +92-^— + 93— M 0 


dqi 


dq<2 


dq 3 


(71) 


The function H will then be convex if M 0 is convex. Since M 0 has been shown to be convex, it is enough 
if H defined by equation (69) satisfies equation (71). For this purpose we observe that 


f = vF f 

dq 1 dq 2 dq 3 


2 ( 7 -l) u 2 5 — 37 
3-7 2 4(7- 1)/? 


F 


(72) 


These equations directly follow from equation (31). Further, equation (71) can be written in the form 

H = IJ{^ + ^ + ^~ F ) dvdI (73) 


which says that the integrand of H defined by equation (73) is also a Legendre transform of F. Using 
equation (72) gives 


dF 3 F OF „ 

q 1 f 92 ^ h 93 a — — F 

dq 1 dq 2 dqs 


In p+ - pu 2 + 20uv - dv 2 - ^ — — 

7—1 3-7 


di- 


5-37 


2(9 -1)J 


(74) 


Equation (58) for In F can be slightly rewritten as 
In F + ^ In d = In p + — -7 + In 


2(7-1) 

Equation (74) then simplifies to 

dF dF dF „ 

Ql "5 b <72 a b 93 ^ — — F 

oq\ dq 2 dq% 


7 - 1 


4(7 ~ 1) 
\A(3 - 9) 


- /?(u - u) 2 - ^-^U/ 


lnF + , 5 3 ? x In/? 


2(9 - 1) 

Hence, the Legendre transform of F is given by 


— F (in 


4 ( 9 - 1 ) 


+ 


E* 


dF 


F = F[\nF + 


5-39 

2(9-1) 


In d — s In 


L v /^(3 — 9 > J 2(9 
4(9-1) 


5-39 ) 
2(9-1)/ 


LaA(3-9)J 2(9 


+ 


!(9- 7 l)}) 


(75) 


Equation (75) is purely a consequence of the Maxwellian distribution and the definitions of < 71 , < 72 , and < 73 . 
A comparison of equations (69) and (70) with equation (75) shows that the //-function and //-flux defined 
by equations (69) and (70), respectively, can be equivalently written as 


*-//*" (e«|£-') 

Hv = f j dvdIv(jT^ t -F^ 


(76) 

(77) 


if we choose 


Ci = - In 


4(9 ~ 1) ] _ 3-9 
l\M 3 - 9 )J 2(9-1) 
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The above analysis reveals a very interesting property of the //-function. First, the //-function is based 
on the Boltzmann //-function and is negative of the density of thermodynamic entropy. Second, because 
of the appearance of In F in its integrand, the //-function is a Legendre transform of another convex 
function. The convexity of H is a consequence of this fact. Further, because of the presence of the 
F\nF term in the //-function, the H - function satisfies the entropy conservation. Thus, the Maxwellian 
distribution is at the root of the convexity of H and the satisfaction of the entropy conservation. 

Finally, we return to the important point mentioned at the beginning of this section, namely, the 
connection between the //-theorem and the entropy condition (eq. (56)). The convexity of the //-function 
implies that 

H (q n+1 ) < H (q n ) 

if the numerical solution <? n+1 at the n + 1 time level is given by 

q n+1 = {l-e)tf) n + 0tf') n (78) 

Here q l and q n could be the values of q at neighboring grid points. Thus, if a numerical method for the 
solution of the Euler equations is of the type used in equation (78), then the //-function will decrease 
with increasing n for that method. 


Extension to Multidimensional Case 

The extension of the results given in the previous section to the multidimensional case is fairly 
straightforward. We will now briefly outline various significant steps involved in extending the above 
analysis to two-dimensional Euler equations: 


dw dgx dgy _ 

dt dx dy 


where 


w = 


p 


pui 


pu 2 

pu 1 

gx = 

P + pu 1 

gy = 

pu 1U 2 



pu 2 


pU\U2 

P + PU2 

. pe . 


_(pe + p)u i. 


,(pe + p)u 2 . 


(79) 


(80) 


As before, equation (79) can be written in the moment form 


T dF dF dF\ n 

*’1^ + V i7r + V 2ir > =0 

at ox ay , 


(81) 


where the moment function 9 is equal to 1, v\, V 2 , and / + /2 j, and F is the Maxwellian. 

Thus, 


F = -p — exp 
I 0 7T 


-P(v 1 - Ui) 2 - 0(V2 - u 2 ) 2 - Y 

I o. 

2- 7 


Io = 


2 ( 7 - 1 )/? 

where I Q is positive for 1 < 7 < 2. Notice that 

w = (*, F) gx = <W1», F) g y = F) 


(82) 

(83) 

(84) 


Equation (84) states that the unknown vector w and the flux vectors g x and g y in equation (79) are the 
moments of the Maxwellian F. We will now show that 


dM 0 , , dM x 

= — — K9x)i = 


dq t 


dqi 


, v _ dMy 
{9y)i ~ 


(* = 1,2, 3, 4) 


(85) 
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<72 = 2u t p 
<73 = 2 u 2 P 
94 = -2/3 


> 



= u\ + ufj 


( 86 ) 


and 


The differentials are given by 


/ oo roo roc 

dv 1 / du 2 / d/F 

/ OO /'OO roo 

dtp / dr> 2 / dl (v\F) 

-oo J —oo JO 

/ OO roo roo [ 

dui / dv 2 / <3/ faF) 

-oo J — oo Jo ) 


dM 0 " 

r 

l 1 ” 

dM x 

= / dtp dv2 dl 
J 

Vl 

dM y 


_^2_ 


dF 


(87) 


(88) 


where, for the sake of notational brevity, only one integration symbol is shown and the limits are also not 
explicitly displayed. Using 


dF J dF , dF dF 

dF ~ ~F P d '’ + dZ[ du ' + &r 2 iU2 + d/3 


d£_F 

dp ~ p 


^=2/3 K-«,)F 


dF 

dU2 


= 2 P(v2 - U2)F 


dF 

dp 


\-{n- ui) 2 - {v 2 - u 2 ) 2 - ^ — — / 
P 2 — 7 


F 


yields 


dF dp 

“^r = h 2p(«i - «i)dui + 2/3 (v 2 - u 2 ) du 2 + 


- («2 - W2) 2 - — / 

2-7 . 


0 “ («i - «l) 2 


dp 


= - 2pui dui -u\dP- 2 Pu 2 du 2 -u\dp + -jj-) 

+ (2/3^i du\ + 2v\u\ dfi) -f (2/3v2 du 2 + 2 V 2 U 2 d/3) 

= d(lnp + 21n/3-/3ii 2 ) + «i d(2/3ui) + v 2 d(2/3u 2 ) 

+ ( /+ y) d(-29) + ^Mu0 


Using the identity 



simplifies the above equation for dF to 


~ =dqi+ vi dq 2 + v 2 dq 3 + + ~ j dq 4 - ^ _ 2 ~ / ^ /d/3 

/ v 2 \ 3 — 2^y 

= dg x + t>i dg 2 - F2 <% + I I + y I d? 4 + ; 2 _ (F, - /) d<7 4 


( 89 ) 


Substituting dF from equation (89) into equation (88) gives 

= J dv\ dv 2 dl 


dM a 
dM x 
dM i 


2 / 


■ 1 " 

- 



.^ 2 . 

- 


/ t> 2 \ 3 _ 2"y 

+ ui d® + v 2 dq 3 + I / + — I dq 4 + — _ ' 


(io ~ /) dq 4 


F (90) 


Again, from observing that 


equation (90) simplifies to 


j (I 0 — I)F dl = 0 


(91) 


dM 0 

r 

’ 1 " 


/ w 2 \ 

d^d x 

— / dv i c/^2 ^7 
J 

^1 


dgi + rq dq 2 + v 2 dq 3 + (/+—] dq 4 

dMy 


_^2_ 


V 2 / J 


(92) 


which implies the validity of equations (85). In terms of the functions M 0 , M x , and M y , the Euler 
equations (eq. (79)) assume the symmetric hyperbolic form 


/ d 2 M 0 dq 3 d 2 M x dq 3 d 2 M y dqj \ _ 

ij \ dq i dq 3 dt d q i dqj dx dqidqj dy J “ 

The proof about the positive definiteness of the matrix P where 

p= a 2 M„ ' 

proceeds along exactly the same lines as before. We first observe that 


(93) 


(94) 


d 2 M, 


°~ I 


dv i dv 2 dl 


dqi + vi dq 2 + v 2 dq 3 + ( / + — j dq 4 


dF 


(95) 


Substituting dF from equation (89) into equation (95) gives 

d 2 M 0 = J dv i dv 2 dl 


dq\ + v\ dq 2 + v 2 dq 3 + I / + — j dq 4 


( v 2 \ 3 _ 2-v 

/ + — I dq 4 + ^ _ (J 0 — /) d(jq 


dq\ + ui d<7 2 


(96) 
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Using / dl (I Q — I) F = 0, equation (96) reduces to two forms: 


*-/■ 


d M 0 = / dv i dl 

4- J dv i dv 2 d/ 


+ v l dq2 + V 2 dqz+ f / + — 1 dq 4 


(u 0 ~ i*) 


and 


«•-/ 


d M 0 — / di>i dv 2 dl dq\ + v\ dq 2 4- t >2 ^93 + 


/ 


4- / dtq dv 2 d/ 


L2-7 

With the use of the easily provable result 


3 7 (^ 2 - //o) (d-74) 2 


7 ~ 1 T v 2 3-7 ^ j 

2^ /+ T + 2^4 7o J 94 

F 


J I 2 Fdldv = j 2I 2 F dl dv 


we obtain 


d 2 M 0 


J dv 1 d ^2 dl |d(?i + i>i d <?2 + t >2 d <73 4- 

- J dv l dv 2 dI [^o(^4) 2 ] F 



dq 4 


(97) 


$M 0 


J dv 1 (*.’2 dl |dgi + V! dq2 + Vi d <?3 + 

+ Y~- J dv i dv 2 dl [lo{du?\ F 


7~ 1 

2-7 



3-27 

2-7 



1 2 


d Q4 


(98) 


If 7 > 1.5, then the positive definiteness of P follows from equation (97); and if 7 < 1 . 5 , then the positive 
property of P follows from equation (98). 

We now come to the derivation of the //-function for the two-dimensional case. As before, the //- 
function is defined as the Legendre transform 


H(F) = H = J (^ ^ 4 / - f) dv i dv 2 dl 

and the fluxes are defined by 


*«=/»> (E«g-*) 


dv i dt >2 dl 


> 

dv\ dv 2 dl 


(99) 


( 100 ) 


From equation (89) it follows that 



Hence. 


1 v — > dF In (3 o 

P 9* ^7 = ln P + [ - 0« + 2/3vi«i + 2/?u 2 u 2 


- 2 /? 


7- 1 


3-27 


i 2 - 1 ,+ Y + 2-7 

ln/3 


/o 


7 - 1 


3-27 


7 

The Legendre transform of F is given by 


= lnp+ -/?(vi - uiY - 0{v 2 - u 2 ) - 2/31 - 

7—1 2-7 7— 1 


ST qi lf-F = F 

^ d Qi 


= F 


ln 


ln 0 


P + - 0{v\ - uif - 0{v 2 - u 2 ) 2 - y— -2(31 

1 -L & '*1 


2-7 


lnp + 


2-7 

- Ul ) 2 - ^2 - W 2) 2 - 7 - - 2 7 


7- 1 


7-1 , ' v ' 1 - 1 ' r ' v ”* I 0 7-1 

In order to express the right-hand side of equation (101) in terms of ln F, we make use of 


In F = ln p + 2 ln (3 + ln 
Combining equations (101) and (102) gives 


2(7 - 1) 


L (2 — 7 ) 71 - J 


(3(v 1 - ui) 2 - 0{v 2 - u 2 f - y 

* O 




2-7 


7 - 1 


2(7 - 1) 


L (2 — y)7rj 


= F\\nF +- — ^ln (3 - C 


where 


7 - 1 


C = + ln 

7- 1 


2(7-1) 


(2 — x)?rj 


The definitions in equations (99) and (100) for H and //-fluxes, respectively, therefore reduce to 


H 

Hyx 


H 


vy J 


= J v\ F^lnF + — — ln/J — Fj dv\ dv 2 dl 


Using equation (104) now shows that 


d H vx OHyy 

dt dx dy 


For this purpose we observe that 


dH dH v 


dt + dx 


+ 


dl L 


f ( dF dF dF\ 

dy = J ll - C + ]nF) \'dr + v 'te +n d^) iv ' iv *dl 

Yl / (I + + ” (Fln/ ’ ),tol *2 dI 


+ 


3- 


(101) 


( 102 ) 


(103) 


(104) 


(105) 


(106) 
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Now, 


[ aj (dF OF dF\ , , Jr 

J l ), {- iH + v ' lh +V2 ai ) dv ' dv i dI 


= 0 


d . d d , ' 

- (p i 0 ) + + - (puM 


\dp d d p d 0 pu\ 80 pu 2 80 

2(7- 1) [dt + dx^ pU ^ + dy^ pU2) 0 dt 0 8 x 0 8 y\ 


7 


7 


2(7-1) 

2-7 


d d d \ 

p Ft ,nll + pui ai [nl, + im2 ry ' nl! ) 

d d d 

— (pin/?) + ^(pu 1 In/?) + — (pu 2 In/?) 


2(7-1) 

= -5^1) / (f + »>| +*^1) 


(107) 


In deriving equation (107) repeated use of the mass balance equation is made. Combining equations (106) 
and (107) gives 


dH dH m dH, 


dt + dx 


+ 


vy 


dF dF\ , , 

+ v\ “ — \-V 2 — dvi dV 2 dl 
dx dy J 


ay =/(‘- c + lnf )(f 


2-7 


■/( 


1 - C + In F - 


6 — 47 
2-7 ' 


dF 


< 9 F 


< 9 F 


01 -TT + U 1 -F-+V 2 -3- ) dvi dv 2 dl (108) 




az 


ay 


Now, equation (102) gives 


6 — 47 

1 + lnF — C — l 0I 

2-7 

— 1 — C + In p + 2 In /? + In 
2(7- 1) 6-47 


2(7-1) 

(2 - ~r)n\ 


0 {v\ - «i ) 2 - /?(v 2 - U2) 2 


/?/ 


2-7 2-7 J 

= 1 — C + In p + 2 In /? + In 

2-7 


2(7-1) 


L(2 -7 )ttJ 


-{ 


1 + In p + 2 In 0 ■ 


7 


■ /?ic 


■}- 


- 0(v 1 - ui) 2 - /3(V2 - «2r - 2 01 
( v 2 \ 

20v\u\ - 20U2V2 - 2/3 I / + — j 


which is equal to the linear combination of collisional invariants 1, tq, V 2 , and I + [v 2 / 2 ). Hence, the 
right-hand side of equation (108) vanishes, thus proving the entropy conservation of equation (105). Once 
again the validity of the entropy conservation has been shown to rest upon the appearance of In F in 
equation (104). The integrand of equation (104) can be expressed as a linear combination of 1, iq, V 2 , 
and / + (f 2 /2). 

From the above analysis the results for the three-dimensional Euler equations 


aw a gx dgy dgz 

dt dx dy dz 


(109) 
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where 


V -f pu\ 
pn\U2 

pU\U% 

.(pe + p)ui. 


pu 3 
pu 1 U 3 

pu 3 W 3 

P + PU3 
.(pe + p)tr 3 . 


/9Wili2 
gy = P + pu 2 
pu 2 W 3 

_(pe + p)u 2 . 


7-1 2 ' 


can be immediately written as 


F dv 1 dv 2 dl 


Then, 


//?\ 3/2 r /i 

( — ) exp —/?(?; — u ) 2 — — 

\nj [ Io\ 


/o U 


J = 5-37 

0 4(7-1)/? 
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The //-function and //-fluxes are defined by 
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and they satisfy the entropy conservation 

dH dH vx dH vv dH vz 

a? + “af + “af + “af = 0 < 116 > 

The convexity of H is a consequence of equation (115) and the convexity of M 0 . 

Review and Discussion 

The Euler equations of gas dynamics are the moments of the Boltzmann equation when the distribution 
function is a Maxwellian. Further, the collision term in the Boltzmann equation vanishes for the 
Maxwellian distribution implying that the Euler equations are the moments of the collisionless Boltzmann 
equation 

dF dF dF dF 

m +v 'te +V2 ^ + V 3 J 7 = 0 
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This is a first-order hyperbolic partial differential equation and can be cast in the strong conservation law 
form by simply taking iq, v 2 , and inside the differentiation symbols. The entire information about the 
Euler equations is compressed in the single equation for the scalar F. For example, any upwind method 
for the collisionless Boltzmann equation becomes an upwind method for the Euler equations by taking # 
moments of the Boltzmann equation. This moment-method strategy has been fully used in reference 2 to 
construct a new class of upwind methods to obtain the numerical solution of the Euler equations. 

The present paper uses the above connection between the Euler and the Boltzmann equations even 
more and shows that the homogeneity of the flux vector, symmetrizability, and the existence and 
construction of the entropy function are all consequences of the Maxwellian distribution. It may be noted 
that the Euler equations can be cast in two forms: the strong conservation law form and the symmetric 
hyperbolic form. The strong conservation law form is obtained when # moments of the collisionless 
Boltzmann equation are taken, and it reflects the physical principle of conservation. The symmetric 
hyperbolic form, which reflects the hyperbolic nature of the Euler equations, is obtained by transforming 
the field vector w to the vector q. At the root of the w-q transformation is the equation 

) d <75 + !^(/ 0 -/)d<75 

When 7 = 5/3, the gas has only translational degrees of freedom and then I and I Q drop out in 
equation ( 111 ) for the Maxwellian distribution F. In such a case the above equation assumes the simple 
form 

dF , , J . v 2 . 

— =dq i + vi dq 2 + v 2 dq 3 + v 3 dq 4 + — dq 5 

The positive definiteness of matrix P (or, equivalently, the convexity of scalar function M 0 ) can be 
proved from the above equation. It is very interesting to observe that the same equation is used to 
establish the convexity of the //-function. It is an interesting property of the Maxwellian distribution 
that the integrand in the definition of the //-function containing In F can be obtained as the Legendre 
transform of F. One of the important results of the present paper is the demonstration of the convexity of 
the //-function based on the Legendre transform and the positive definiteness of the matrix P. Therefore, 
the Maxwellian distribution F and the transformed variables play a fundamental role in the theory 
of Euler equations. Just as the field vector w naturally arises when the Euler equations are cast in the 
strong conservation law form, the variables q naturally arise when these equations are transformed to the 
symmetric hyperbolic form. 

The //-function, which is a slight modification of the Boltzmann //-function, is the negative of specific 
entropy of thermodynamics. Because of the convexity of the //-function, it is possible to design a method 
satisfying the entropy condition, that is, a method for which total H in the computational domain 
decreases with time. A decrease in total H therefore corresponds to the existence of entropy-producing 
mechanisms in the numerical method. It should therefore be possible to obtain nonoscillatory solutions 
to the Euler equations by controlling the decrease in H. This fact, together with the observation that 
a decrease in H is physically meaningful, in that it avoids expansion shocks and is easily extended to 
multidimensional flows, gives credence to the view that it would be more profitable to construct total 
//-diminishing (THD) methods instead of the well-known total variation diminishing (TVD) methods. 

Finally, the ability of the moment-method strategy to tackle a system of equations by taking moments 
of a single scalar equation can be used in various ways. First, for a system of conservation equations 
it is enough to develop a TVD-lgke criterion for a single scalar equation having a physical meaningful 
connection with that system. Second, in this framework upwind methods for the three-dimensional flows 
can be constructed without the usual splitting of a three-dimensional problem into three one-dimensional 
problems. 
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